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INTRODUCTION 


An emerging trend in the analytical design of aircraft is the integration of all appropriate 
disciplines in the design process (refs. 1,2). This means not only including limitations on the 
behavior of the design from the various disciplines, but also defining and accounting for 
interactions so that the disciplines influence design decisions simultaneously rather than 
sequentially. The integrated approach has the potential to produce a better product as well as a 
better, more systematic design practice. In rotorcraft design (the rotor in particular), the 
appropriate disciplines include aerodynamics, dynamics, structures, and acoustics. The purpose 
of this paper is to describe a plan for developing a helicopter rotor design optimization procedure 
which includes the above disciplines in an integrated manner. 

Rotorcraft design is an ideal application for integrated multidisciplinary optimization. There are 
strong interactions among the four disciplines cited previously; indeed, certain design parameters 
influence all four disciplines. For example, rotor blade tip speed influences dynamics through 
the inertial and air loadings, structures by the centrifugal loadings, acoustics by local Mach 
number and air loadings, and aerodynamics through dynamic pressure and Mach number. All of 
these considerations are accounted for in current design practice. However, the process is 
sequential, not simultaneous, and often involves correcting a design late in the design schedule. 

Applications of rigorous and systematic analytical design procedures to rotorcraft have been 
increasing, especially in the past five years. Procedures have accounted for dynamics (refs.3-9), 
aerodynamics (refs. 10-11), and structures (ref. 12). Generally, these applications have only 
considered single-discipline requirements, although in reference 5, dynamic and structural 
requirements were considered together, and in reference 6, dynamics and aeroelastic stability 
were combined. Integrated multidisciplinary applications to rotorcraft are largely nonexistent. 

In early 1985, several occurrences led to an excellent opportunity at the NASA Langley Research 
Center to address the multidisciplinary design problem for rotorcraft. The Interdisciplinary 
Research Office was established and charged with the development of integrated 
multidisciplinary optimization methods. Nearly concurrently, the Army Aerostructures 
Directorate at Langley established the goal of improving rotorcraft design methodology by 
"discipline integration." Close cooperation between the NASA and Army organizations led to 
initial plans for a comprehensive, integrated analytical design capability. By 1986, a group of 
NASA/Army researchers had formed a committee and began detailed planning for this activity. 
The committee, designated IRASC (Integrated Rotorcraft Analysis Steering Committee), has 
now completed the bulk of the planning and has formulated the approach described in this paper. 


I 


The development of an integrated multidisciplinary design methodology for rotorcraft is a three- 
phased approach. In Phase 1, the disciplines of blade dynamics, blade aerodynamics, and blade 
structures will be closely coupled, while acoustics and airframe dynamics will be decoupled from 
the first three and will be accounted for by effective constraints on the other disciplines. In Phase 
2, acoustics will be integrated with the first three disciplines. Finally, in Phase 3, airframe 
dynamics will be fully integrated with the other four disciplines. This paper deals only with the 
Phase 1 approach and includes: details of the optimization formulation, design variables, 

constraints, and objective function as well as details of discipline interactions, analysis methods, 
and methods for validating the procedure. 
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LANGLEY PARTICIPANTS IN THE ACTIVITY 


The work described in this paper represents the combined efforts of a team of researchers and 
managers at the NASA Langley Research Center and Army Aerostructures Directorate involved 
in analysis and optimization of rotorcraft. Shown on the left side of figure 1 are the members of 
the IRASC (Integrated Rotorcraft Analysis Steering Committee). Shown on the right side of the 
figure are the technical contributors. This team includes experts in the areas of rotorcraft 
aerodynamics, dynamics and aeroelasticity, structures, and acoustics, as well as optimization and 
sensitivity analysis. 
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Figure 1 
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OBJECTIVE OF THE EFFORT 


Figure 2 represents the Charter of the IRASC activity. The unique features of the work are the 
emphasis on integrating the disciplines and explicitly accounting for the interactions among 
disciplines. While the team does not intend to develop new or improved analyses in any of the 
included disciplines, the latest developments in rotorcraft disciplinary analysis will be used. 
Finally, it is a goal of the team to stimulate activity in the rotorcraft (and the aircraft) community 
in the general area of multidisciplinary design integration and the use of formal Mathematical 
Programming in design work. 


Develop and validate an integrated multidisciplinary 
design capability involving aerodynamics, dynamics, 
structures, and acoustics which leads to improved 
design practices and improved rotorcraft. 


Figure 2 
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FUNDAMENTAL DECISION/STRATEGY 


The rotorcraft optimization effort at Langley has a three-phased plan (fig. 3). The initial phase, 
which is well underway, will decompose the rotor design problem. This will be accomplished by 
driving the design through the integration of the disciplines of aerodynamics, dynamics, and 
structures while satisfying additional constraints imposed by airframe dynamics and rotor 
acoustics. The latter constraints will account for the influences of the airframe response and 
acoustics behavior on the overall optimization process. Phase 2 of the design plan will include 
rotor acoustics as a discipline inside the optimization loop; that is, integrated with aerodynamics, 
structures, and dynamics. Finally, in Phase 3 of the effort, all five disciplines will be fully 
integrated inside of an iterative design loop. This paper focuses primarily on the Phase 1 
activity. 


Phase 1 

• Decompose the design problem 

• Blade aerodynamics, dynamics, structures 

• Airframe analysis/modeling and acoustics 

• Develop representation of acoustics and airframe 
influences on blade aerodynamic, dynamic, 
structural optimization 

Phase 2 

• Bring acoustics into optimization loop 
Phase 3 

• Bring airframe dynamics into optimization loop 


Figure 3 
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FOCUS NO. 1 - WHITE PAPER 


Publication and critique of a "white paper” (fig. 4) detailing the plan for rotorcraft optimization at 
Langley is viewed as a prime focal point of the activity. The paper, which is presently in draft 
form, will present the goals of the design effort, as well as the approach and validation plan. The 
approach will discuss strategy, analytical methods, and discipline couplings. The validation 
procedure will document test problems to be used to examine the fidelity of the specific 
discipline tools used, as well as the overall optimization procedure. Finally, the white paper will 
be disseminated to industry for critique and a workshop is planned to consolidate feedback to the 
plan. 


• Goals of the activity 

• Approach and plan 

• Formulation of optimization strategy 

• Governing mathematics (varying detail) 

• Definition of interactions 

• Analysis methods (codes) to be used 

• Validation methods - test problems 

• Status 

• First draft written 

• To be critiqued by industry 


Figure 4 
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FOCUS NO. 2 - ROTOR DESIGN 


A practical rotor design will be the end result of the Langley optimization effort (fig. 5). The 
purpose of this "fidelity check" includes not only an overall test of the methodology, but also a 
measure of each discipline’s modeling effectiveness. The test problem for Phase 1 will contain a 
rotor task, mission scenario, and challenging design requirements. Simulation models describing 
aerodynamic, structural, and dynamic systems will be formulated to allow for key 
interdisciplinary couplings. Design variables, constraints, and a specific objective function will 
be identified. As the design activity progresses, so will the verification of individual discipline 
models. Finally, the optimum design will be wind-tunnel tested for fidelity of the entire Phase 1 
process. 


• Apply Phase 1 method to rotor design 

• Define test problem 

• Rotor task and mission 

• Design requirements 

• Generate math models 

• Aerodynamic 

• Dynamic 

• Structure 

• Incorporate couplings 

• Formulate optimization problem 

• Objective function 

• Design variables 

• Constraints 

• Validate design methodology 

• Discipline tools 

• Overall 


Figure 5 
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INTEGRATED ROTORCRAFT OPTIMIZATION DEVELOPMENT PLAN 


Figure 6 depicts the general sequence of tasks that will lead to a fully integrated rotor blade 
aerodynamic/dynamic/structural optimization procedure which also accounts for acoustic and 
airframe dynamic influences. The dynamic optimization work is building on the work described 
in references 5, 6, and 9. The rotor aerodynamics activity has been separated into two parts. The 
first is aerodynamic performance optimization which is a continuation of the work described in 
reference 10. The second is an integration of aerodynamic loads analysis with dynamics - a 
procedure wherein the local airloads can be adjusted by varying the planform dimensions and 
twist of the blade to reduce dynamic response. A merger of the rotor performance optimization 
with the airload/dynamics optimization will yield a fully integrated aerodynamic/dynamic 
procedure. The rotor structural optimization is a continuation of the work of reference 12. A 
merger of all the aforementioned procedures, with the acoustic and airframe constraint interfaces, 
will lead to the fully integrated Phase 1 procedure. The resulting capability will be applied to a 
rotor test article to validate the procedures. 
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DEFINITION OF THE OPTIMIZATION PROBLEM 


As indicated in figure 7, the next section of the paper consists of details of the integrated 
rotorcraft optimization problem. Included are descriptions of the following: the objective 

function (the quantity to be minimized for obtaining an optimum design); the design variables 
(dimensions and other parameters of the design); constraints (a set of behavioral or characteristic 
limitations required to assure acceptable and safe performance); and definitions of the 
interactions among the disciplines. 


• Objective function 

• Design variables 

• Constraints 

• Interactions among disciplines 


Figure 7 
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CONTRIBUTORS TO OBJECTIVE FUNCTION 


In formulating the objective function, a number of indicators (fig. 8) were considered. Basically, 
the indicators fell into two categories - penalty types such as weight and cost; and performance 
types such as vibratory loads (shears and moments); and required horsepower. It was decided, 
for the purpose of the Phase 1 work, to choose performance type quantities for the objective 
function; specifically, transmitted vertical vibratory hub shear and horsepower required at several 
flight conditions. Blade mass will be prevented from becoming excessive by enforcing an upper 
limit constraint. Cost is not explicitly accounted for in the formulation. 


• Blade mass 

• Cost 

• Vibratory loads 

• Required horsepower 


Figure 8 
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FORM OF OBJECTIVE FUNCTION 


The objective function is a linear combination of the main rotor horsepower at five flight 

conditions plus the transmitted vertical vibratory hub shear at a frequency of N times the blade 

angular speed (where N is the number of blades). As shown in figure 9, the objective function 

contains weighting factors K through K which will be assigned based on a mission criterion 

1 6 

to be determined. The five flight conditions referred to above include hover, cruise, high speed, 
maneuver, and maximum range. The speed and load factor specifications for these conditions 
are given in the lower portion of the figure. 


• Linear combination of main rotor horsepower at 
five flight conditions and transmitted vertical 
hub shear 

•F = Kj HF^+ K 2 Hf£+ K 3 HP 3 + K 4 Hfy+ K 5 HP 5 + KgS 


Flight 

condition 

Description 

Velocity 

(Kts) 

Load 

Factor 

1 

Hover 

0 

1.0 

2 

Cruise 

140 

1.0 

3 

High Speed 

200 

1.0 

4 

Maneuver 

120 

3.5 

5 

Max. Range 

95 

1.0 


Figure 9 
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BLADE MODEL AND DESIGN VARIABLES 


Figure 10 is a depiction of the rotor blade model to be used in the Phase 1 optimization activity. 
Also shown in figure 10 are the design variables which are defined in figure 11. The blade model 
can be tapered in both chord and depth. The depth can be linearly tapered from root to tip. The 
chord is constant from the root to a spanwise location (referred to as the point of taper initiation) 
and is linearly tapered thereafter to the tip. Design variables which characterize the overall shape 
of the blade include the blade radius, point of taper initiation, taper ratios for chord and depth, the 
root chord, the blade depth at the root, the flap hinge offset, and the blade maximum twist. 
Tuning masses located along the blade span are characterized by the magnitude and locations. 
Design variables which characterize the spar box beam cross-section include the wall thicknesses 

at each spanwise segment and the ply thicknesses at 0° and +45.° Additional design variables 
include the number of rotor blades, the rotor angular speed, and the distribution of airfoils. 



Figure 10 
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SUMMARY OF DESIGN VARIABLES 


Description 

Symbol 

Tuning mass at location i 

Mj 

Spanwise location of i-th mass 

Xi 

Wing box dimensions 

ti , t 2 , t3 

Ply thicknesses 

*45 » *o 

Depth of blade at root 

h r 

Ratio of blade depths at tip and root 

= ^r^t 

Maximum pre-twist of blade 

xmax 

Percent blade span where taper begins 

r 

Width of blade at root 

^r 

Airfoil distribution 


Hinge offset 

e 

Blade angular velocity 

Cl 

Number of blades on rotor 

N 

Blade radius 

R 

Ratio of root chord to tip chord 

! = 0 p/Cj 


Figure 11 
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CONSTRAINTS OVERVIEW 


As previously described, the Phase 1 activity is based on integrating the blade aerodynamic, 
dynamic, and structural analyses within the optimization procedure. The acoustics and airframe 
dynamics analyses are decoupled from the first three disciplines and their influences are 
expressed in terms of constraints. Accordingly, the total set of constraints is made up of two 
subsets as indicated in figure 12. The first subset consists of constraints which are evaluated 
directly from the first three disciplinary analyses and are a direct measure of the degree of 
acceptability of the aerodynamic, dynamic, and structural behavior. The second subset 
represents indirect measures of the satisfaction of constraints on the acoustics behavior and the 
requirement of avoiding excessive vibratory excitation of the airframe by the rotor. 


• Aerodynamic 

• Dynamic 


• Structure 


Evaluated directly 
from analysis 


• Acoustic 

• Airframe 


Expressed indirectly 
► as constraints on 

disciplines in top group 


Figure 12 


14 



SUMMARY OF CONSTRAINTS 


The constraints are summarized in figures 13 and 14. The first two constraints are for 
aerodynamic performance and require that for all flight conditions, main rotor horsepower not 
exceed available horsepower and that airfoil section stall not occur at any azimuthal location. 
The next nine constraints address blade dynamics. The first requires that the blade natural 
frequencies be bounded to avoid approaching any multiples of rotor speed. The next five impose 
upper limits on the blade vertical and inplane loads, transmitted hub shear, and hub pitching and 
rolling moments. The next three dynamic constraints are an upper limit on blade response 
amplitude, a lower limit on blade autorotational inertia, and finally, the aeroelastic stability 
requirement. The structural constraints consist of upper limits on box beam stresses, blade static 
deflection, and blade twist deformation. The acoustic constraints are expressed as an upper 
bound on the tip Mach number and tip thickness to limit thickness noise; and an upper bound on 
the gradient of the lift distribution to limit blade vortex interaction (BVI) and loading noise. The 
effective airframe constraints are expressed first as a separation of the fundamental blade inplane 
natural frequency in the fixed system from the fundamental pitching and rolling frequency of the 
fuselage. Second is a bounding of the blade passage frequency to avoid the proximity to any 
fuselage frequency. The final constraint is an upper limit on the blade mass which will avoid any 
designs which satisfy the constraints at the expense of large mass increases. 


Constraint Description 

Form of Constraint 

Comments 

Main rotor horsepower 
Airfoil section stall 

HR HP avail for 
i-th condition 

Cp ^Cpmeix 

For 5 flight 
conditions 
Enformed at 
12 azimuthal 
locations 

Blade frequencies 
Blade vertical load 
Blade inplane load 
Transmitted in-plane 
hub shear 

Hub pitching moment 
Hub rolling moment 
Blade response amp. 
Autorotational inertia 
Aeroelastic stability 
Wing box stresses 

fil*fj*fiu 
Vjk<Vmax 
Hjk'SHmax 
Xk< X max 
Yk< Ymax 
P k < Pmax 
Rk*sRmax 
qk'SQmax 
Imjrj 2 > a 
Re(A,)<-e 
R<1R = Tsai -Hill 
criterion 



Figure 13 
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SUMMARY OF CONSTRAINTS - CONCLUDED 


Constraint Description 

Form of Constraint 

Comments 

Blade tip deflection 
Blade twist 

W'SW max 
6< Omax 


Blade tip Mach no. 
Blade thickness 

M<Mmax 
h< hmax 

Limits 

thickness 

noise 

Blade lift distribution 

dC|/dx<Smax 

Limits BVI 
& loading 
noise 

Ground resonance 
Rotor/Airframe 
frequency coupling 

|ft-W M |«W af 

f|-SN£2«f u 

Effective 

airframe 

constraint 

Blade mass 

M'SM upper 



Figure 14 
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ANALYSIS ASPECTS 


The analytical tools (summarized in fig. 15) must provide technical fidelity in phenomena 
prediction, as well as connectivity between disciplines. The areas of aerodynamics, dynamics, 
and structures will utilize codes to predict response, as well as sensitivity information. The 
constraint-providing disciplines of acoustics and airframe dynamics have the analysis task of 
defining the impact of the design on acoustic energy and fuselage response. 

The aerodynamic analysis for rotor performance prediction will include a hover momentum/strip 
theory code for hover and climb applications (ref. 13). The CAMRAD analysis will be used for 
forward flight and maneuver performance. In order to assure that the latest developments in 
inflow analyses are available, some modularity will be provided in the inflow modeling based on 
recent fidelity assessments (ref. 14). 

Rotor dynamics will utilize CAMRAD for forced response calculations. Finite element modeling 
(ref. 15) and modal analysis (ref. 16) will form the tools for the dynamic tuning before the global 
analysis predicts the final blade loads, response, and rotor stability. 

The structural codes involve a combination of beam analysis and laminate analysis. The beam 
analysis (e.g., ref. 12) is applied to the blade planform model. The laminate analysis will be 
applied to one or more cross-section models. The beam model consists of equivalent stiffness 
and masses from which displacements and forces are computed. The internal blade structure is 
represented by cross-section models to calculate resultant stresses associated with each beam 
model segment. The laminate analysis then uses these stresses to determine critical structural 
margins of safety. 

The effectiveness of imposing Phase 1 acoustic constraints will be quantified by using the 
WOPWOP code (ref. 17), with appropriate loading inputs from CAMRAD. Low frequency 
loading, thickness, and BVI noise will be generated from this analysis. 

Airframe dynamics constraints for Phases 1 and 2 will result from fixed-system frequency 
predictions and will neglect hub motion. Phase 3 of the effort will involve finite element 
modeling and impedance tailoring to effect favorable rotor-body coupling in the design process. 
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ANALYSIS ASPECTS 

Hover momentum/strip theory 
Global forward flight code 
Inflow modules 
Finite element modeling 
Modal code 

Global analysis/forced response 
Eigen analysis for stability 
Beam models 
Laminate model 
WOPWOP code to verity 
acoustic acceptability 
Fixed system frequency prediction 
Finite element models - NASTRAN 


i 

i 

Figure 15 


Aerodynamic 

Dynamic 

Structure 

Acoustic 

Airframe Dynamics 
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INTERDISCIPLINARY COUPLINGS 


Phase 1 of the Langley rotorcraft optimization effort will utilize several design variables which 
have historically been significant drivers of disciplinary phenomena. In addition, other variables 
are being included to provide other unexplored design opportunities. Figure 16 shows our first 
attempt to quantify the interactions among the disciplines through the design variables. For 
example, rotor tip speed has driven past rotor designs based solely on acoustics, performance, or 
dynamics. This variable also influences blade structural integrity and fixed system response to 
transmitted loads. This provides the strong interdisciplinary coupling for tip speed shown in 
figure 16. There are variables, such as blade twist, which can strongly influence some disciplines, 
such as aerodynamics, while not perturbing others (e.g., structures) and other variables such as 
hinge offsets which, heretofore, have not greatly influenced conventional rotor design. 

It is a significant part of the current design methodology effort to explore not only the obvious 
strong design variable couplings, but also to address those variables which may provide design 
synergism for multidisciplinary design goals. This may provide a design key for missions which 
have not been accomplished with today’s rotorcraft. 


Variable 

Acoustics 

Aerodyn. 
(Perf & Loads) 

Dynamics 

Structures 

Fuselage 

Dynamics 

Airfoil Dist. 

S 

S 

W 

w 

W 

Planform 

S 

S 

S 

S 

S/W 

Twist 

W 

S 

S 

w 

W 

Tip speed 

S 

S 

S 

s 

S 

Blade number 

S 

W 

S 

w 

S 

Stiffness 

W 

S 

s 

s 

S/W 

Mass dist. 

W 

w 

s 

s 

S/W 

Hinge offset 

W 

w 

s/w 

w 

S/W 


S = Strong interaction 
W = Weak interaction 


Figure 16 
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INTEGRATED AERODYNAMIC/DYNAMIC/STRUCTURAL OPTIMIZATION 

OF ROTOR BLADES 


Figure 17 is a flow chart which explains how the integrated procedure will function. The current 
set of design variables (summarized in fig. 11) will be input to design variable preprocessors 
which will generate input tailored for each analysis. For example, box beam cross-sectional 
dimensions will be used to compute values of El and GJ for use in the dynamic analysis. Each 
disciplinary analysis will be carried out using the preprocessed data along with the necessary 
input which is the product of other disciplinary analyses. For example, as the flow chart shows, 
the dynamic response requires airloads from the aerodynamic loads analysis. The appropriate 
output from each disciplinary analysis is collected in a module which computes the objective 
function (fig. 9) and constraints (figs. 13 and 14). The next major step is the sensitivity analysis 
to calculate derivatives of the objective function and constraints with respect to the design 
variables. At this stage, all the information is available for the optimizer module (based on the 
program CONMIN (ref. 18)) wherein the values of the design variables are updated. The above 
steps are repeated until a converged design is obtained. Convergence requires that the objective 
function is minimized and all constraints are satisfied. 


Current design variables 
and fixed parameters 


Design variable preprocessors 


Aerodynamic 
analysis 
(loads) 


Airloads. 


Dynamic 

analysis 


Response 

stability 


Dynamic 


loads 


Structural 

analysis 


4 


Aerodynamic 

analysis ->|Objective function and constraints [< 
I (performance) T 


Stresses 

deflections 


HP req'd 
trim 


I 


Sensitivity analysis 

zzzuzzz: 


Optimizer 


Updated design variables 


Figure 17 
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RESULTS OBTAINED TO DATE 


As indicated in figure 18, progress has been made in the areas of aerodynamic performance 
optimization, dynamic optimization, optimum placement of tuning masses for vibration 
reduction, and structural optimization. Selected results from these activities arc highlighted in 
the next portion of the paper. 


• Aerodynamic performance optimization 

• Dynamic optimization 

• Optimum placement of tuning masses 

• Structural optimization 


Figure 18 
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RESULTS - AERODYNAMIC PERFORMANCE OPTIMIZATION 


A Mathematical Programming technique has been developed to minimize the hover horsepower 
for a helicopter with a specified design gross weight operating at a specified altitude and 
temperature (fig. 19). A conventional design approach is a two-step iterative method. The first 
step is design for optimum hover performance by varying taper ratio, point of taper initiation, and 
twist until the rotor blade configuration with the lowest hover horsepower is obtained. In the 
second step, this best hover design is modified by changing the root chord to meet forward flight 
and maneuverability requirements. The Mathematical Programming approach uses the same 
performance analyses as the conventional approach, but couples a general-purpose optimization 
program to the analyses. The conventional and Mathematical Programming approaches have 
been used to define the blade configuration which provides the lowest hover horsepower and 
satisfies forward flight and maneuverability requirements. The figure also summarizes results for 
the final design variable values and the main rotor horsepower required for hover from each 
approach. The Mathematical Programming approach produced a design with more twist, a point 
of taper initiation further outboard, and a smaller blade root chord than the conventional 
approach. The Mathematical Programming design required 25 less hover horsepower than the 
conventional design. Most significantly, the Mathematical Programming approach obtained 
results more than ten times faster than the conventional approach. 


Conventional 


h i 

" ' 

r 

Hover 

Forward 

flight 



Mathematical 

programing 


CDesigner^) 


Hover 


Forward 

flight 


Optimizer 


Objective 


Design 
variables 


Design time 


Hover horespower 

1558 hp 

1533 hp 

j 

[ Twist, deg 

-12 

-15 

J 

1 Percent taper 

.80 

.91 

1 

l Taper ratio 

3.0 

3.1 

1 

{ Root chord, ft 

2.3 

1.78 

ne 

5 weeks 

2 days 


Figure 19 
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RESULTS - DYNAMIC OPTIMIZATION 


A rotor dynamic optimization problem is summarized in figure 20. Blade mass is the objective 
function. Upper and lower bound constraints are placed on the first five natural frequencies 
(elastic modes only) of the blades to separate them from the excitation frequencies. Also, a 
lower bound constraint is imposed on the blade autorotational inertia. A stress constraint is used 
to guard against structural failure due to blade centrifugal stress, and side constraints are imposed 
on the design variables to avoid impractical designs. The design variables are the box beam spar 
wall thicknesses, the magnitudes of tuning masses, the blade taper ratio, and the box beam height 
at the blade root. The program CAMRAD and CONMIN are used for the blade modal analysis 
and the optimization, respectively. A sensitivity analysis calculates analytical derivatives of the 
objective function, autorotational inertia and stress constraints, along with finite difference 
derivatives of the frequency constraints. 


Objective function - blade mass W 

W = Wjj + W Q Wfo = box beam mass 

W Q = nonstructural mass 

Constraints 

Frequency: f j( < fj<fj u fj ►first three 

Autorotational inertia: lead-lag 

^ \ first two flap 

X W.r.^ ? a a: minimuim rotary 
Vj = 1 J I j inertia 

Stress constraints: o/centrifugal < aa 
Design variables - Box beam cross-sectional dimensions 

Tuning masses 
Blade taper ratio 


Figure 20 
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RESULTS - DYNAMIC OPTIMIZATION (CONCLUDED) 

Optimum designs for minimum mass rotor blades have been obtained for both rectangular and 
tapered blades. The optimum designs are compared with an existing baseline blade denoted the 
’reference’ blade in figure 21. The reference blade is based on an actual flight article and is 
described in more detail in references 5 and 9. The figure also shows the box beam wall 
thickness distributions for the rectangular blade using 30 design variables (t^ , t^, t^ at ten spanwise 

locations), and for the tapered blade with 42 design variables (ten lumped masses, and h^ and 

are the additional design variables). Blade mass reductions of four to six percent (compared to a 
baseline or reference blade) have been achieved without violating the imposed constraints. The 
optimization process tends to shift the wall thickness distribution outboard for both designs due 
to the presence of the autorotational inertia constraint. The tapered design requires more 
outboard mass shift, but this is easily accomplished within the stress constraints. 


Design variables: 


X.h = 7 X. 



OPTIMUM HORIZONTAL WALL 
THICKNESS DISTRIBUTIONS 


.01 2i 
.0 

Horizontal 
wall 

thickness^) oo6| 

fl .004 
.0021 



Reference 
Optimum (rect) 
Optimum (tapered) 


l 


r-r 




.83 4.4 8.8 13.2 17.6 

Radial distance, r 


22 . 


Blade weight 
(lbs) 

Percent reduction 


Reference y 


98.3 


Optimum 


Rectangular 
93.6 

4.7 


T apere d | 
92.1 

6.2 


Figure 21 
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RESULTS - OPTIMUM PLACEMENT OF TUNING MASSES 


The design problem (shown in figure 22) is to find the best combination of tuning masses and 
their locations to minimize blade root vertical shear without a large mass penalty. The objective 
function is a combination of vertical shear and the sum of the tuning masses. Constraints arc 
placed on the frequencies to avoid resonance. The strategy employed reduces the total shear as a 
function of time during a revolution of the blade. 

The example problem (figure 23) is a beam representation of an articulated rotor blade. The 
beam is 193 inches long with a hinged end condition and is modeled by ten finite elements of 
equal length. The model contains both structural mass and lumped (non-structural) masses. 
Three lumped masses are to be placed along the length of the beam. The strategy was applied to 
a test case of two modes responding to three harmonics of airload. The figure compares the 
blade vertical root shear s(t) plotted versus azimuth for the initial and final designs. The peaks 
on the initial curve have been reduced dramatically. For example, the oscillatory (1/2 p-p) load 
was reduced from 78 lbs to 0.6 lbs - nearly two orders of magnitude. 


cb « 



• Design goal - Find optimum combination of 
masses and their locations to reduce blade 
root vertical shear 

• Method - Formulate optimization procedure 

• Use masses and locations as design variables 

• Minimize - 

• Blade root vertical shear 

• Added mass 


Figure 22 
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RESULTS - OPTIMUM PLACEMENT OF TUNING MASSES (CONCLUDED) 


cb “ Blade 
I model 


M. M, 
1 £. 


M 


Results based on two 
vibration modes and 
3,4,5, per rev airloads 


Initial design 
Final design 



Initial 

Final 

X 1 

(in) 

135.0 

153.0 

X 2 

(in) 

154.0 

154.0 

X 3 

(in) 

174.0 

154.0 

m; 

(Ibm) 

5.2 

7.9 

Mo 

(Ibm) 

6.5 

9.7 

m 3 
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Figure 23 
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RESULTS - STRUCTURAL OPTIMIZATION 


A blade structural optimization procedure (fig. 24) applicable to metal and composite blades has 
been developed in which the objective function is blade mass with constraints on stresses in the 
spars and in the skin, twist deformation, and autorotational inertia. The design variables are the 

total spar thickness and for the composite blade the percentage of ±45° plies (the remaining plies 

assumed to be at 0°). This procedure is described in detail in reference 12, and additional 
applications of the methods are also given in reference 12. 


• Objective function: Blade mass 

• Constraints: Stress in skin and spars, 

twist deformation, 
autorotational inertia 

• Design variables: Spar thicknesses, 

% of ±45°plies 



Figure 24 
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RESULTS - STRUCTURAL OPTIMIZATION (CONCLUDED) 

Two examples rotor blade design (fig. 25) were generated using the structural design 
methodology. Both designs use the UH-60A Black Hawk titanium spar blade as a baseline. 
First, a titanium spar blade design was generated. A maximum twist deformation of 3.1 degrees 
was selected. The structural constraint requires that the calculated stresses not exceed the 
allowable material strength which is assessed using the Tsai-Hill failure criterion (ref. 12). The 
autorotational capability is satisfied by requiring the mass moment of inertia to be at least 19000 
2 

in-lbs-s per blade. 

As shown in figure 25, the minimum metal spar thickness which satisfies all of the constraints is 
.130 inches and corresponds to a blade weight of 207 pounds. The actual UH-60A titanium spar 
is .135 inches thick and weighs 210 pounds. The new titanium spar design is only three pounds 
different from the actual UH-60A blade, demonstrating that the mechanics of the design 
methodology can produce blade designs similar to conventional design processes for the same 
design allowances and material. A second design was developed using a single T300/5208 
graphite/epoxy D-spar. The composite design also satisfied the required constraints and the 
minimum weight design had a .105 inch thick spar with 20 percent of the plies oriented at +45 
degrees. The composite blade weighed 163 pounds which represents a 21.5 percent reduction 
over the actual UH-60A blade. 


Twist constraint 



Autorotation constraints 

A| = Jjfiu > 17 

1100 hp 



Material strength constraints 

\l/ 



R = Tsai-Hill 
failure criterion 


Results 


Parameter 
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Metal Comp, 
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Ti 

Ti 

Gr/E 

Spar thick., in. 

0.135 

0.130 

0.105 

45° plys, % 

— 

... 

20 

Weight, lbs. 

210 

207 

163 

Margin (1-R) 

— 

0.103 

0.000 

Twist, deg. 

— 

0.92 

2.55 


% of 45° 
plys (all 
others Or) 



Design variables 

t (thickness 
of spar) 


Figure 25 
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VALIDATION OF PROCEDURES 


The process of validating the optimization methodology (fig. 26) involves substantially more 
than evaluating the success of the final design. Specifically, the analyses used in optimizing the 
rotor during the Phase 1 effort will be examined for predictive fidelity and design technique 
validation. The usefulness of the basic tools involves not only accuracy of analysis, but also a 
reliable parametric sensitivity capability. Several opportunities are currently available to assess 
the fidelity of the analyses. For example, rotor performance, dynamics, and acoustics predictions 
need accurate inflow distributions for various flight conditions. Recent experimental efforts 
(e.g., ref. 19) and code validations (ref. 14) are helping to provide confidence in the available 
inflow models. Rotor geometric design variable sensitivity (e.g., effect of taper on performance) 
which was reasonably well-known for past rotor designs, is being re-examined in light of recent 
correlation anomalies for high-speed flight. Acoustic source mechanisms and modeling validity 
are also being examined (ref. 20), especially for parametric sensitivity of the acoustic energy to 
rotor state. Structural coupling mechanics are being exploited in new rotor designs to assess the 
structural tailoring benefits while satisfying structural integrity requirements (ref. 22). 

Proof of the fidelity of design techniques is crucial to the overall design optimization effort. For 
example, aerodynamics and dynamics interact so strongly in rotor design that basic aeroelastic 
tailoring efforts must be validated. Such a validation effort is being undertaken at Langley, as 
well as other research centers (ref. 22). Also, since rotor speed is a strong driver for aeroelastic 
response, a program to assess variable RPM designs is underway at Langley. The object of this 
effort is to define the benefits and limitations of an aerodynamically and dynamically designed 
multi-speed rotor. In addition to design techniques, which capitalize on the strong effects of 
certain design variables, small variances in other blade characteristics may impede the practical 
operation of even conventional designs. Hence, the ability to accurately predict even these 
secondary phenomena is important for the design effort. For example, track-and-balance 
sensitivity experiments and studies are being undertaken which can lead to a practical design 
capability to eliminate blade-to-blade variability effects. 
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VALIDATION OF PROCEDURES 

• Basic predictive tools 

• Rotor inflow fidelity 

• Acoustic source mechanisims 

• BVI 

• Loading 

• Thickness 

• Composite structional couplings 

• Parametric taper initiation studies 

• Design technique validations 

• Blade dynamic tailoring 

• Variable RPM designs 

• Blade-to-blade variability effects 


Figure 26 
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OVERALL DESIGN VALIDATION 


For the overall Phase 1 validation effort (fig. 27), the Langley team has chosen a rotor task which 
requires maneuverability, speed, and efficiency. Specifically, the rotor mission must be 
accomplished with minimum power and vibration while satisfying predefined acoustic, stability, 
and fuselage dynamics requirements. 

The assessment of the Phase 1 design methods will involve model rotor hover and wind tunnel 
tests. The models (a baseline and an advanced design) will be aerodynamically and dynamically 
scaled. Provisions for varying key design parameters are necessary to complete the validation 
process. In other words, the tests need to quantify not only a minima, but the gradients around it. 
The testing possibilities include a series of 1/5-scale model rotors, mounted on a variable drive 
system and tested in hover and simulated forward flight in a tunnel which can eliminate many 
testing "excuses" such as inappropriate Reynolds, Mach, and Froude Numbers. The Langley 
TDT is the candidate facility for the major segments of the validation process. 


• Definition of mission profile at gw/altitude/power 

• Maneuverability 

• Speed 

• Efficiency 

• Minimize objective function while satisfying constraints 

• Wind tunnel model rotors 

• Baseline 

• Advanced design with some parametric variability 

• Aerodynamicaly and dynamically scaled 


Figure 27 
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CONCLUDING REMARKS 


This paper has described a joint NASA/Army initiative at the Langley Research Center to 
develop optimization procedures aimed at improving the rotor blade design process by 
integrating appropriate disciplines and accounting for important interactions among the 
disciplines. The activity is being guided by a Steering Committee made up of key NASA and 
Army researchers and managers. The committee, which has been named IRASC (Integrated 
Rotorcraft Analysis Steering Committee), has defined two principal foci for the activity: a 
"white paper" which sets forth the goals and plans of the effort; and a rotor design project which 
will validate the basic constituents, as well as the overall design methodology for 
multidisciplinary optimization. 

This paper has described the optimization formulation in terms of the objective function, design 
variables, and constraints. Additionally, some of the analysis aspects were discussed and an 
initial attempt at defining the interdisciplinary couplings was described. At this writing, some 
significant progress has been made - principally in the areas of single discipline optimization. 
Results were given in the paper which represent accomplishments in rotor aerodynamic 
performance optimization for minimum hover horsepower, rotor dynamic optimization for 
vibration reduction, and rotor structural optimization for minimum weight. 


• Defined objectives and procedures for integrated 
rotorcraft design 

• Formalized Army/NASA team - 
Researchers and Managers 

• Goals - Two foci 

• White paper to be critiqued by industry 

• Rotor design to validate procedures and tools 

• Disciplinary optimization and coupling 
methodologies already yielding useful results 

• Rotor aerodynamic performance 

• Rotor dynamics 

• Rotor structure 


Figure 28 
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